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Abstract 



An asymmetric generalization of the zero-temperature Glauber model 

on a lattice is introduced. The dynamics of the particle-density and spe- 

►^ \ cially the large-time behavior of the system is studied. It is shown that the 

system exhibits two kinds of phase transition, a static one and a dynamic 
one. 



1 Introduction 

In recent years, reaction-diffusion systems Irave been studied by many peo- 
ple. As mean-field techniques, generally, do not give correct results for low- 
dimensional systems, people are motivated to study exactly-solvable stochastic 
models in low dimensions. Moreover, solving one- dimensional systems should in 
principle be easier. Exact results for some models on a one-dimensional lattice 
have been obtained, for example in [1-12]. Different methods have been used 
to study these models, including analytical and asymptotic methods, mean field 
methods, and large-scale numerical methods. 

Some interesting problems in non-equilibrium systems are non-equilibrium 
phase transitions described by phenomenological rate equations, and the way the 
system relaxes to its steady state. Kinetic generalizations of the Ising model, for 
example the Glauber model or the Kawasaki model, are such phenomenological 
models and have been studied extensively [13-18]. Combination of the Glauber 
and the Kawasaki dynamics has been also considered [19-21]. 

In this article, we want to study an asymmetric generalization of the zero- 
temperature Glauber model on a lattice with boundaries. There are also sources 
(or sinks) of particles at the end points of the lattice. We study the dynamics 
of the particle density, and specially the large time behavior of the system. In 
the thermodynamic limit, the system shows two kinds of phase transition. One 
of these is a static phase transition, the other a dynamic one. The static phase 
transition is controlled by the reaction rates, and is a discontinuous change 
of the behavior of the derivative of the stationary particle density at the end 
points, with respect to the reaction rates. The dynamic phase transition is con- 
trolled by the injection- and extraction- rates of the particles at the end points, 
and is a discontinuous change of the relaxation time towards the stationary 
configuration. 

2 Asymmetric Glauber model at zero tempera- 
ture 

In the ordinary Glauber model, the interaction is between three neighboring 
sites. Spin flip brings the system to equilibrium with a heat bath at temperature 
T. A spin is flipped with the rate /i := 1 — tanh -^ if the spin of both of 
its neighboring sites arc the same as itself, and is flipped with the rate A := 
1 -I- tanh -^ if the spin of both of its neighboring sites are opposite to it. At 
domain boundaries, the spins are flipped with unit rate. So the interactions can 
be written as, 

AAA -> A^HA and 000^0^0 ^ 
A9A->AAA and 0^0^000 A 
AA0^A00 and 00A^0AA 1 



where spin up and spin down are denoted by A and 0. One can interpret an up 
spin as a particle, and a down spin as a hole. At zero temperature, the Glauber 
dynamics is effectively a two-site interaction | |2^ : 

where all the above processes occur with the same rate. 

One can consider the following interactions, as an asymmetric generalization 
of the zero-temperature Glauber model. 



AA 



u 



^'-[m : (^' 



AA 



V 



If M ^ w, the above system has left-right asymmetry. The above system on 
an infinite lattice has been investigated in p3[ , where its n-point functions, its 
equilibrium states, and its relaxation towards these states are studied. It can 
be easily shown that the time evolution equation for the average densities of the 
system with the above interactions are the same as that of a system with the 
following interactions, where diffusion is also present: 

A 

u-A (5) 

i; - A 

A' 




- A' . (6) 



V 



u-A' 



Consider a lattice with L sites and an asymmetric zero-temperature Glauber 
dynamics as the interaction. The rates of injection and extraction of particles in 
the first site (final site) are a and a' (b and &'), respectively. The time evolution 
equations for the average densities are then 

(hk) = -{u + v){nk) +u{nk^i) +v{nk+i) for ky^l,L 

{ill) = a— {a + a' + v){ni) + v{n2) 

(til) = b-{b + b' + u){nL)+u(nL-i). (7) 

First, let us calculate the profile of average densities at large times. At large 
times, the system goes to its stationary state, and the time-derivatives of the 
left-hand sides vanish. One can see then that the solution to the above system 
is 

(nfc(oo))=a + /3(^)'. (8) 



Putting this in M), one can easily find a and (3. In the thermodynamic hmit 
{L -^ C!o), the solution becomes 

ha' - ah' /'u\'^-L 

u> V 



a ba — ab /u\ 

^"fc(°°^> = ^T^+ (a + a')(h + h' + u-v) \v) 



a + a' {a + a'){h + h' + u-v) 

b ab' — a'h /u\*^^-i 

(nkioo)) — r; + T, tttt ; 7 ~ i u < v. (9) 

^^ ^' 6 + 6' (6 + 6')(a + a +«-") ^^^^ 

It is seen that the density profile is flat at the left end {k <^ L) for u > v and 
its value is independent of the reaction rates. But as v exceeds u, the density 
profile acquires a finite slope, proportional to ln(u/v). This is the static phase 
transition previously mentioned. 

Now return to the dynamics of the system. Let us write the homogeneous 
part of (^ as 

{hk)^hi{ni), (10) 

and find the eigenvalues and eigenvectors of the operator h. One finds 

Exk = —{u + v)xk + uxk-i + vxk+i, k^l,L 

Exi = —{a + a' + v)xi + VX2, 

Exl = —{b + b' + u)xL+UXL-l, (11) 

where the eigenvalue and eigenvector have been denoted with E and x, respec- 
tively. The solution to these equations is 

Xk^az'l + f3z^, (12) 

where Zi^s satisfy 

u 
E = -{u + v)+vz+-, (13) 

and 

v{azf + Pz^)- {E + a + a' + v) {azi + f3z2) = 
uiaz^-^ +(3z^~'^)^{E + b + b' + u){az^ + (3z2) = 0. (14) 

To have nonzero solutions for x, these last two equations should be dependent, 
the criterion for which is 

(u + ziSa) [vz^^^ + z^Sb) - {u + Z2Sa) ivz^+^ + z^5b) = 0, (15) 

where ( |l3[) has been used to eliminate E, and Sa ;= a + a' — u and Sb := b + b' ~v. 
Defining 

Zi '■= Zi\ — 
V u 

Sa 

A := 



Sb , , 

B := ^, (16) 



(Esh is simplified to 

^-(i+i)(l _,_ AZ){1 + BZ) - Z^+^{1 + AZ-'^){1 + BZ-'^) = 0. (17) 

The eigenvalue E satisfies 

E^-{u + v) + ^/mi{Z + Z~^). (18) 

Two obvious solutions of the equation (|l7| ) are Z = ±1. But these generally 
don't correspond to eigenvalues and eigenvectors. In fact for these solutions, Z 
and Z~^ are the same, so that (n3) should be modified to 



Xfe = (a + /3fc)(±l)* 



(19) 



and it is not difficult to sec that these do not fulfill the boundary conditions 
unless a — P ^ 0. Equation (|l^) can be written as a polynomial equation of 
order 2(L + 1), and hence has 2L more roots in addition to Z = ±1. For these 
2L roots, if Z is a root Z~^ is another root, and these two correspond to one 
eigenvalue and one eigenvector. So the L x L matrix h does have L eigenvalues 
and eigenvectors. 

For A = B = 0, (hTu is very simple and its nontrivial solutions are 






1 < s < i. 



(20) 



and their inverses. All of these are phases and the real-part of the corresponding 
eigenvalues satisfy 



Re(£') < -{u + v) + 2V™ cos 



L + 1 



<0. 



(21) 



The maximum of the real-part of the eigenvalues determines the relaxation time 
toward the stationary average-density profile. That is 



(22) 



1 



u + V — 2vuucos 

In the limit i — > oo, this is simplified to 

r = {u + V — 2\/wJ)^^ . 
The general solution to (|^) is seen to be 

2 ._ /u\(fe-™)/2 



(23) 



[nk 



W) = E rTT^''''f^^^^*^^"'^°^^ (iJ) 



sm 



TTsk \ . / nsm 
sin 



L + 1 



L + 1 



(24) 



Now consider the general case. The equation (^^ can be written as 

G{Z) := F{Z) - F{Z-^) = 0, (25) 



where 

F{Z) -.^ Z-^^+^\l + AZ){l + BZ). (26) 



If Z is a phase, it satisfies (17[) provided F{Z) is real. Consider the phase of 
F{Z) for unimodular Z. We have 

cj)[F{z)] = -{L + l)(t){Z) + 0(1 + AZ) + 0(1 + BZ), (27) 

where denotes the phase of its argument. As the phase of Z is changed from 
to 27r, the change of the phase of i^(z) is 

A0[F(Z)] = -2tt{L + 1) + A(/)(l + AZ) + A(/)(l + BZ). (28) 

For |A|, |i?| < 1, the phase changes of (1 + AZ) and 0(1 + BZ) is zero, as Z 
moves on the whole unit circle. So, 

A0[F(Z)] = -27r(L + 1), for|A|, \B\ < 1. (29) 

But this means that the phase of F{Z) will be an integer multiple of it for at 
least 2{L + 1) points on the unit circle. So all of the solutions of ( p^ ) are still 
phases, although they may be not uniformly spaced on the unit circle. Two of 
these are Z = ±1. The remaining 2L points correspond to L eigenvalues and 
eigenvectors for h. One concludes that for \A\, \B\ < 1, 

T={u + v-2y/^cose)^'^, (30) 

for some 9. Specially, at L — ^ cxd, the relaxation time is the same as the relax- 
ation time for A = B = 0, that is the same as (p^). 

If, for example, |^| > 1, then the total phase change of (1 + AZ) is no 
longer zero. It is 2tt. One may then lose two of the roots of the unit circle. 
Note that the mere fact that A0[F(2')] = —2ttL does not mean that there are 
just 2L solutions of { pj\) on the unit circle, since the phase of F(Z) needn't be 
monotonic. To find values of A and B, for them the number of the solutions of 
( [l7| ) on the unit circle is 2L or 2(L— 1), consider the function G{Z) at the points 
Z — ±1. Increasing \A\ or \B\, two of the roots on the unit circle tend to 1 or 
— 1, and then move out of the unit circle and on the real line. At the point that 
this occurs, either G"(l) or G"(— 1) become zero, as there will be multiple roots 
at ±1. So the criterion for each change (losing 2 roots of the unit circle) is that 
either G"(l) or G"(— 1) become zero. The curves in the AB plane, corresponding 
to these changes are 

L{l + A){l + B) + 1-AB = 0, G'(1) = 0, (31) 

and 

L{l^A){l-B) + l-AB = 0, G'(-1) = 0. (32) 

These curves divide the plane into six regions: 
I: all of the solutions are phases. 



II: 2L phase solutions, 2 real negative solutions. 

Ill: 2(L — 1) phase solutions, 4 real negative solutions. 

IV: 2L phase solutions, 2 real positive solutions. 

V: 2(i — 1) phase solutions, 4 real positive solutions. 

VI: 2(i — 1) phase solutions, 2 real negative solutions, 2 real positive solutions. 

In the above, by real solutions it is meant real solutions besides the trivial 
solutions ±1. Note, however, that not all of this plane is physical. The physical 
region is that part of the plane, which corresponds to nonnegative values for 
the injection and extraction rates. Returning to the definitions of A and B, 
equation (Q), it is seen that 

V 

: <-) 

The point {Aaiin, -Smin) itself is on the curve 

^111111-0111111 = 1- ("54) 

One can see that, unless u — v, part of the physical region is in the region 
IV, where two of the solutions of (|l^) are real positive (and of course inverse 
of each other). This makes the maximum of the real part of E larger than 
— {u + v) + 2y/uv, and correspondingly the relaxation time larger than (123). 

In the thermodynamic limit L ^ oo, the regions of the AB plane are greatly 
simplified. In fact the curves corresponding to G'(±l) = become 

A = -l, or B = -l, G'(1) = 0, (35) 

and 

A=l, or -8 = 1, G"(-1) = 0. (36) 

In the case either |yl| or \B\ are greater than 1, the real roots of ( [l7|) arc simply 

Z = -A,-A-\ \A\>1 (37) 

and 

Z^~B,~B~\ \B\>1 (38) 

Then, if for example A is negative and less than —1 and A < B, the maximum 
real part ofE'is ~{u + v) — y/uv{A + A^^), and the relaxation time of the system 
is 

T = [u + V + V^{A + A-^)]-\ (39) 

which is greater than (E3h. This is a phase transition which occurs at A = — 1. 
For A > —1, the relaxation time is constant, (ESh. For A < —1, it is A- 
dependent. The minimum of A is — -y/u/w, for which one of the eigenvalues of h 



become zero, and the relaxation time becomes infinite. The same effect is seen 
for B < —1. As mentioned before, for u ~ v (the ordinary zero-temperature 
Glauber model) no part of the region IV is in the physical region, and this 
transition does not occur. The phase transition discussed here, is the dynamical 
phase transition mentioned before. 
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